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Abstract 

We develop a molecular model of non-uniform deformations within the framework of liquid 

crystal rubber elasticity. We show that, similarly to the uniform case, the theory is not sensitive to 

the molecular details of polymer liquid crystal involved and the resulting elastic free energy is quite 

universal. The result of this work is the general expression for the free energy of deformations, which 

combines the effects of large non-symmetric amne strains in the rubbery network and gradients of 

curvature deformations of the director field, F ~ A T {Vn} 2 A. We derive the molecular expressions 
Ph. 

for the 5 principal independent elastic constants governing non uniform deformations in the presence 
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of elastic strains. They also depend on the polymer step length anisotropy and - most strikingly - 
have an overall negative sense. We therefore predict that in some circumstances, especially at large 
elastic deformations A, these new terms may overpower the usual, positive Frank elastic moduli of 
^ ■ the underlying nematic structure, as well as the coupling in nematic elastomers of uniform relative 

1 rotations of the director and the elastic matrix. In this event highly distorted polydomain textures 
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n(r) would be favoured. 



PACS: 61.40K - 61.30B - 46.30C 



I. INTRODUCTION 



A large amount of experimental and theoretical work has been invested in recent years into studies of liquid 
crystalline elastomers and gels. Elastomers (as opposed to hard and brittle resins) are weakly crosslinked, percolating 
networks of polymers chains which retain a significant molecular mobility of their strands. Rubber, therefore, is a 
material with a very low shear modulus /i (about 10~ 4 — 10~ 5 times that of conventional elastic solids) and thus 
deforms at practically constant volume. Its weakness is matched by its capacity to withstand very large deformations, 
commonly 100's % and more, thus allowing large amounts of elastic energy to be stored in the material. 

Nematic liquid crystals (and liquid crystal polymers) also represent an intermediate state between conventional 
liquids and solids. They flow under applied stress, but they have a long range orientational order and thereby a 
curvature elasticity associated with deformations of this order. When these two remarkable physical systems are 
combined in the same material, i.e. a weakly crosslinked network is formed of mesogenic polymer chains capable of 
spontaneous orientational ordering, a qualitatively new state of matter emerges. Elastomers, being nominally solids, 
display a high molecular mobility since connected crosslinks are distant from each other and the chains of mesogenic 
monomers have much freedom. The axis of orientational symmetry breaking (the nematic director n) is mobile too 
and responds to imposed elastic strains. We have, therefore, a uniaxial solid where the preferred direction can be 
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altered when the body is elastically deformed; n becomes an independent elastic variable. Such systems resemble the 
so-called Cosserat medium, where the internal torques are allowed and, therefore, elastic stress can be non-symmetric. 

For most deformations n is anchored in the elastic medium by the crosslinking. Its rotation requires the input of 
some external work, but at large strains it may jump in direction. For some deformations the director is (theoretically) 
freely mobile, although still intimately coupled to the network, and deformations occur without resistance, that is 
without applied stress. Thus nematic elastomers seem to exhibit a qualitatively different elasticity from other solids 
and, in the non-linear regime, display nematic and mechanical transitions hitherto unknown. Elastic non-linearity is 
an important motivation for developing a molecular theory of nematic elastomers: rubbery networks are capable of 
extremely large extensions and continuum elastic models are bound to break down well before elastomers do. 

Another important motivation is that a specific picture of how the mesogenic chains spontaneously change their 
equilibrium shape on entering the nematic phase is unnecessary for a detailed description of most liquid crystal 
elastomer properties. This allows a universal description of nematic elastomers, a modest extension of the classical 
molecular theory of conventional rubbers, independent of molecular models of liquid crystal polymers where there is 
much less agreement. 

Experimental studies of liquid crystal elastomers began more than a decade ago and various side-chain |^-||] 
and main-chain |6|-^) materials have been investigated. However, it was quickly discovered that if the mesogenic 
polymer is crosslinked in the isotropic phase, the nematic phase obtaining on cooling down through the clearing 
point T„i inevitably has a scattering, highly nonuniform equilibrium texture of the director. Although the question 
of quenched random disorder created by misoriented crosslinks has some fundamental interest, for most practical 
effects and applications a uniform monodomain birefringent elastomer was preferrable. Several methods of obtaining 
such systems have been developed, notably crosslinking in a magnetic field || and two-step crosslinking with stress 
applied in the intermediate state pC| . Experiments then showed the effects of coupling of the mobile anisotropy axis 
to the elastic strain field in such monodomain nematic elastomers, for example the switching of the director by the 
perpendicular extension has been observed [pd| . Interested reader can find more information in the recent review 
articles HQ. 

There has been a substantial effort in developing both the continuum (applicable only for small deformations) and 



the full molecular theory of nematic elastomers, which are also reviewed in 13 1. Much of this effort has been successful 
in describing the existing and predicting new physical effects. However, all these theoretical models addressed only the 
problem of uniform deformations and there has been no adequate theory describing the effect of director curvature. At 
the same time, there are many mechanisms of director deformation, for example due to surface anchoring, disclinations, 
or in transient regimes. Therefore, it is important to know how director curvature is affected by elastic deformations of 
the rubber. The present paper is devoted to this particular point: after a brief overview of basic concepts and results 
of the classical theory in the next Section, we derive the full non-uniform nematic rubber elastic free energy, which 
depends on the products of elastic strains tensors, chain step length tensors and on the second power of Vn (only in 
chiral - cholesteric - systems can we expect linear gradients of n ) . In the last Section we examine the general free 
energy form via simple examples and compare its effect with that of standard Frank elasticity. The main conclusion 
of this paper is that some of the new non-uniform rubber elastic terms are negative and, in some circumstances, may 
favour equilibrium non- uniform textures n(r). 
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II. BASIC THEORY OF UNIFORM NEMATIC ELASTOMERS 



The first linear continuum picture of the local anchoring of the director with respect to the elastic matrix is due to 
de Gennes |l5|| , the energy density being (phenomenologically) 

E rot = ifoi [(H - u>) x n] 2 + b 2 n ■ e • [(fi - w) x n] . (1) 

The local rotation of the elastic medium is £1, given by the antisymmetric part of the infinitesimal deformation tensor, 
e^: Qi = tijk&jk- is the rotation of the director n about an axis parallel to u>, the change in the director being for 
small rotations 8n = [uj x n]. The first term in Eq.([l]) expresses the penalty for the relative director-matrix rotations. 
It is quite unlike usual ncmatic Frank elasticity, which depends on gradients of director rotation. In ncmatic solids 
uniform rotations are also penalised and the degeneracy of the local direction of orientational order is removed. The 
62 term is the first hint that it is not only rotations of the anchoring matrix that can rotate the director, but that the 
symmetric component of shear strain is coupled to the director as well. 



More recently, molecular models of nematic elastomers have been constructed [[16 17 1, which led to an understanding 
of mechanical critical points, memory of crosslinking, shifts in phase equilibria and stress-strain relations. The essential 
anisotropy of polymer chains leads to a straightforward modification of the conventional Gaussian rubber elasticity 
theory. The theory can be further developed, as in conventional networks, to account for junction point fluctuations 
jl7j and other effects, but to understand the startling new effects visible in nematic elastomers these elaborations 
are not necessary. The theory is based on a single parameter, the anisotropy of the polymer strand shape, i.e. on 
the ratio of the mesogenic chain persistence lengths along and perpendicular to the nematic director, £\\/£±. One 
should note that this ratio strongly depends on the molecular nature of the polymer (£||/£j_ ~ 1.5 in side-chain liquid 
crystal polymers |jl|[l!|, while in a main-chain system one can find t\\/l±_ ~ 10 or more [Q) and that different values 
are predicted by the different theoretical approaches (which include freely-jointed chains, worm-like persistent chains, 
rotationally-isomeric chains, etc.). The happy situation in the theory of liquid crystal elastomers is that this single 
parameter i\\/(± is directly related to the macroscopic sample shape and can be easily measured simply by changing 
temperature through the nematic-isotropic transition || and observing the spontaneous sample shape change. After 
such a measurement the theory does not have any free parameters and should make accurate predictions. 

The basic theory assumes that in a nematic monodomain with director n°, a given polymer chain span between 
connected crosslinking points, at the moment of crosslinking, is R° . If this chain is long enough (so that we obtain 
an elastomer, not a resin), this end-to-end distance has a Gaussian distribution: 

F (R°) ~ Det[/°.]- 1/2 ex P (-^(l^R^j (2) 

(Summation over repeated indices has been assumed). The matrix £®j of effective step lengths defines the chain shape 
parallel and perpendicular to the director n° for a uniaxial phase, that is 



<i?°i?°) = -U% , (3) 



where C is the chain contour length. The effective step lengths of the random walk, and thus the overall average 
shape, are functions of the nematic order parameter. For a uniaxial nematic one can write: 

hi = Z±Sij + - £±) rii rij , (4) 
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where the explicit functional form (\\(Q) and £±(Q) is dependent on the particular model of the liquid crystal polymer 
chain. We shall not need this information to describe the rubber-elastic effects. 

Next comes the affine deformation assumption, an assumption that pervades network theory. Junction point 
fluctuations are damped by connectivity and in the case of very high crosslink functionality compel crosslinking 
points to move exactly in geometric proportion to the whole sample. We use the simplifying but unnecessary affine 
deformation assumption and define the current network span to be Ri = X%jRj with Ay the macroscopic (Cauchy) 
strain tensor of the whole block of rubber. The polymer strand end-to-end vector probability is given by a distribution, 
-P(R), differing from Po(R°) in Eq.(||) in that the (^ij) -1 tensor, describing the current chain shape, depends on the 
current state of the nematic order after the deformation has taken place. Taking the usual quenched average for 
each network strand F e i = — (fc^T In P(R))p (ro) , one obtains the elastic free energy density describing uniform 
deformations of liquid crystal elastomers [glj : 

F el = ^N x k B T Tr[£° -X T -t- 1 -X] , (5) 

where N x is the number of crosslinks per unit volume (see [ p~7[ for insignificant corrections due to the junction point 
fluctuations) and we have used = Eq.(|J). There is also an additional term which arises from the 

normalization of -P(R) (see fl7|Jj2|j ). This term depends on the magnitude of the nematic order parameter Q, which 
could, in principle, also be affected by elastic strains. At sufficiently low temperatures, away from the clearing point 
T n i, the degree of polymer chains anisotropy £\\/£± is strongly constrained by thermodynamic forces and it is quite safe 
to assume Q = const and consider only the rotations of the nematic director n in response to elastic strains in Eq. (^|) , 
so that £ij would be just a rotated version of the matrix £^ (see the review article [|l3| for a detailed discussion). 

The elastic free energy (|j) showns a rich behaviour, based on the fact that the elastic strain tensor enters this 
expression in a generally non-symmetric form, unlike in any other solid material (or isotropic rubber), where one 
always has F e i ~ {A T ■ A }. Antisymmetric components of strain couple to the director rotation away from no (given 
by £°) to n (the principal axis of the current step length tensor £). In particular, after implementing the limit of 
infinitesimal deformations A = 1 + e , one obtains molecular expressions for de Gennes' coupling terms in Eq.(Q): 

h = N x k B T L- ; b 2 = N x k B T \— (6) 

(a positive coupling b 2 corresponds to prolate polymers, those with mesogenic units in or aligned parallel to the 
backbone and thus with £« > £±). Several new physical effects have been predicted with the help of Eq.(^), notably 
discontinuous nematic transitions driven by an imposed elastic strain (also observed experimentally Jll| ), reorientation 
by external electric or magnetic fields and the so-called 'soft elasticity'. However, the above arguments are applicable 
only in the case of uniform director rotations and uniform elastic strains. Applications of Eq.(||), however appealing 
and powerful they might look, are severely hindered by its neglect of curvature deformations which are quite common, 
even in elastomers, due to surface anchoring, disclinations and domain walls in equilibrium, or during orientational 
transitions. 

It is the purpose of this paper to derive the non-uniform elastic free energy, by a modification of Eq.(^|), that 
describes the coupling between curvature deformations and elastic strains in nematic elastomers. 
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III. NON-UNIFORM DIRECTOR FIELD 



We shall assume, as a simplifying starting point, that the director distribution before deformation, no, which is 
implicitly present in the initial chain step length tensor A, is uniform. In other words, the polymer network has been 
formed in a uniform monodomain nematic state, or brought to its present state in a sufficiently strong external field. 
Another simplifying assumption will be that the present temperature is sufficiently below the clearing point T n i, so 
that the magnitude of the nematic order (and the related chain step length anisotropy £\\/£±) is not changed during 
the deformation and only director rotation takes place. Both these assumptions are not crucial for the arguments 
below and a straightforward generalization to account to these effects is possible. However, these two factors, Vno 
and £ij(Q), bring additional degrees of freedom and geometric constraints, which make the resulting elastic free energy 
very difficult to read. We prefer to present the main line of derivation in a more clear, albeit slightly less general form 
and, therefore, consider no = const and £\\/£± — const. 

In order to determine the elastic response of a random polymer network to a non-uniform deformation field one has 
to find the statistical weight P(R) of configurations of a given polymer strand in the presence of such a field. This 
statistical weight for a chain with a fixed contour length, C, is determined by the general tensor of chain persistent 
lengths £ij, see Eqs.(^), (||). Then the quenched averaging should be performed, exactly as in Eq.(|J), with respect to 
the initial state P(R°), which is characterised by the corresponding initial £^ and which we assumed to be uniform. 
Therefore, we need to describe the configurations of a liquid crystal polymer chain in the presence of director curvature 
deformations £ij(n, Vn). After expansion in powers of small gradients (long- wavelength limit) we should then have 
in a uniaxial non-chiral nematic: 

e v ~ 4? + K ijabcd{n)V a n b V c n d + ... , (7) 

which tells us that the shape of the polymer coil is slightly altered when the mesogenic units of this chain are subjected 
to a spatially non-uniform nematic mean field potential represents the uniaxial step length tensor of the uniform 
system). The statistical weight of such a chain with its end-to-end vector R affected by the affine elastic deformation 
of the network is, as before, 

P(R)~e X p^-A { X.^}^r.i(n,Vn){A • R°}^j . (8) 

Performing the quenched averaging with the (uniform) initial probability distribution P(R°), Eq.(||), one obtains 

F e i = ^N x k B T Tr[£° -A T -l' 1 -A] 

-^ B TTr[f°.A r -{...VnVn...}.A] , (9) 

where is the uniform liquid crystal elastomer elastic energy Eq.(|^). Note the minus sign at the non- uniform part 
of F e i, which is due to the inversion £~ x from Eq.(^|) and which will appear to represent a genuinely negative free 
energy contribution. 

How can one find £ij(n, Vn)? Obviously a specific model of the liquid crystal polymer chain must be employed 
and the result could be very different for main-chain mesogenic polymers (best modelled by a persistent worm-like 
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chain) and for side-chain materials (their backbone configuration is adequately described by the much simpler freely- 
jointed rod model). Since the majority of existing liquid crystal elastomers are made of side-chain polymers, we shall 
concentrate on this case. We present a calculation of the average end-to-end distance (RiRj) of the uniaxial nematic 
polymer in the non-uniformly distorted director field n(r), within a freely-jointed chain model. Later we shall also 
briefly outline the similar calculation for the main-chain work-like nematic polymer. 

Let us choose the starting point of the chain trajectory as the origin of coordinate system, then the position of the 
a th monomer on the chain is given by = a^]" =1 where a is the physical length of the monomer and u, y is 

the tangent vector of the monomer number v on this chain. The direct product of end-to-end vectors is, therefore, 
RR = a 2 EL'=i uMu{i/} 

, where N — C/a is the number of monomers on the chain. For the freely-jointed chain 
there is no correlation between orientations of different monomers, and ) for v ^ v', and the shape of the 
polymer coil is determined by 

JV 

(R i R j )=a 2 J2(4 V} ^ } ) (10) 
i/=i 



(we take Latin indices to represent the vector components, i,j = 1, 2, 3, while Greek indices number monomers along 
the chain, a, v = 1, ... , JV. For the freely-jointed chain there is a trivial relation between the monomer length, step 
lengths and the backbone order parameter Q: a = ^(£\\ + 2£±) ; Q — — £±)/a. In the second expression in Eq. (]To|) 
we explicitly expose the fact that the orientation of the principal axis of the uniaxial average (uu) depends on the 
position of the corresponding monomer, numbered v. This is the crux of the problem since itself depends on all 
the preceeding monomers orientations and hence on the nematic director n(r) on all these locations. The formerly 
simple problem of a freely-jointed random walk has become a higher order Markov process, the v th step depending 
on all the previous steps. This problem can be solved systematically in powers of (Vn) 2 . 

The end-to-end vector R of a chain with N monomers can also be written in a recurrent form, R = N ~ 1 } + au^ N ^ 
and thus the average required in Eq. (|l0|) is 

(R i R j ) = (ry-%?-») + a*(ul N ^) (11) 

(cross terms (r^^^u^^) vanish in a non-chiral system for obvious symmetry reasons). The last term in ( |Tl| ) 
is equal to ^a(Sij£± + [£\\ — £±]rii(r^ N ~ 1 ^)nj(r^ N ^ 1 ^)^ where it is explicitly noted that the directions of the N th 
link are biased according to the director r^r^ -1 }) at its starting point r^ -1 ^. Recognizing that in the long- 
wavelength limit the difference between n(r^) and n(r'f A '^) is small (which implies the upper cut-off for the curvature 
wave vectors, |q| <C |R| _1 , the inverted network mesh size), we perform the gradient expansion: rii 
n*(r {1} ) + (r^-^ • V)n l (r< 1 >) + i^-^-^r^- 1 } : VVjii.frf 1 !) + ... , where one should treat n^rW) as the local 



current director n. Substituting this into Eq.flllD we obtain 

{RiRj) = (4 N ~ 1} rf- 1} ) + \a{l^ + ft - i ± ]n inj ) (12) 

+i<zft| - e±]{r^ lf ~ 1} r} N ~ 1} ) (V fc n<)(Vinj-) + -n^JuW^ + J^V^n* 



2 2 



(we have eliminated cross terms like (rii(r^ N 1 ^ • Vn^)) which vanish for a non-chiral system). So far Eq.([L2h is exact 
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at O(Vn) 2 and the problem of correlations has been set back to (rj N ^ ^ r^ N ^). This equation can now be iterated to 



G 



) in terms of (ri 



,{N-2} W-2} 



} and gradients of n, and so on. Each step of such iteration generates 



successively higher powers of [in — £±](Vn) 2 . Since there are no linear terms, ignoring all powers greater than (Vn) 2 
is equivalent to taking the uniform-n value for {r { k v} r\ v} ) [see Eq.@ with N = v\. In this way we obtain 
3 



hi = -^(RiRj) = e±S i:j + [£\\ - ix\n%nj 



(13) 
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-aN\£u 



- £_l)(ii • Vn,)(n • Vn ) + e±(V k n i )(V k n j ) 



where the extra power of N has appeared due to the summation X^=i v = ^N(N + 1) of all 
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) terms in the 



iterated Eq.([l2]). Inverting the matrix (|l^) is trivial because we should only keep the lowest (second in this case) 
powers of Vn; in this limit of small gradients we obtain 



r. 1 



rnk 



(14) 



\u\ 

where is the local "uniform" step length tensor depending on n, and the expression in square brackets should be 
directly taken from Eq.([l3|). 

All that remains is to substitute this inverse of the non- uniform step length tensor in the general local rubber-elastic 
free energy density Eq.(^|) and integrate by parts in order to convert the result to the form consistent with {Vn} 2 
only [instead of having the second-derivative terms, suggested by the (|l3|)]. This integration by parts is quite tedious 
and strictly speaking, since the elastic strain tensor is in general a function of coordinates, will generate derivatives 
of A too. However, it is common in elasticity theories to neglect the gradients of strains (which correspond to second 
derivatives of displacement in infinitesimal models). We, accordingly, shall ignore the terms of the form {VA 2 Vn} 
in favour of the coupling to the uniform part of strains {A 2 (Vn) 2 } , which bears a superficial similarity with the 
nematic Frank free energy and rubber-elastic energy (|^) . In this case the algebra is straightforward and we obtain the 
main result of this paper, the elastic free energy density of non-uniform director deformations in nematic elastomers 
[compare with Eq.(^)]: 

1 



AF Kl = 



-«i Tr 



— —Kn Tr 

2 

— -Ka Tr 

2 d 



K4 Tr 



Tr 



£° ■ \ T ■ {(n ■ V)m (n-VK) A 
£°- X T - {VfcTuVknj} - A 
e° ■ \ T ■ {mrij} ■ X ] [(n-V)n] 2 
1° ■ X T ■ {mnj} ■ X (V k ni)(V k ni) 
1 • A T ■ {m(n ■ V)n/cVfcrij + Uj(n • V)n k \7 k m} ■ X 



(15) 



-Tr 



'°-A T 



{ni(n ■ V)rij + nj(n ■ V)rii} • A 



div n 



where the elastic constants have the form 



«i = -pk B Ta- 
6 



Vi 



K-3 



-pk B Ta- 
6 



M 3 



^2 = -zpk B T a- 
6 



k 4 = -pk B Ta- 
6 
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(16) 
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K 5 = -pk B Ta ^ , 

and where p = N X N is the total number density of monomers in the system. 

There can be different representations of these elastic constants, using chain step lengths, physical monomer length, 
or the backbone order parameter. Experimentally, perhaps the most easily accessible parameters of the material are 
the monomer size a (simply from its chemical structure) and the aspect ratio of the nematic polymer chain A = t\\/t±_ 
(from neutron scattering or from the sample shape change on isotropization), which is also the single parameter of 
the "uniform" theory pl[ and determines all elastic instabilities. In these variables, the five non-uniform nematic 
rubber-elastic constants take the form: 

1 a (A-l) 3 1 a (A-l) 2 

Kl = 7;P k sTa A , A —r ■ k 2 = ^pk B Ta ^ — ^ ; (17) 

1 . rp 2 U-lr 1 a (A-l) 2 

«3 = —pkB-L a ; K4 = —pkBi a 



(A- 


1)3 


A(A- 


f2) 


(A- 


1)3 


A 2 (A 


+ 2) 


(A- 


I) 2 


A(A- 


f2) 



2 r A 2 (A + 2) ' * 2 r D A 2 (A + 2) 

k 5 = ^pk B Ta' 

In order to complete these equations, an additional factor of £± has been pulled out of the tensor t in the traces in 
Eq. (|l5|) , so that the initial step length tensor there should now be regarded as dimensionless, & — Sij +(A— 1) n°n°. 
For freely jointed rod polymers £± = 3a/ (A + 2). One can immediately notice that the dimensionality of the elastic 
constants n is energy per length, the same as the Frank constants of a nematic. 

The reader should be reminded that the above derivation is performed explicitly within the framework of freely- 
jointed rod model for the polymer chain. In the most general case there are two more non-uniform nematic rubber 
elastic terms, arising from the persistent correlations along the chain. 



-KfiTr 

2 6 



£° ■ \ T ■ {Sij} ■ X] [(n-V)n] 2 - ]-K 7 Jr I u • A 1 ■ {%} • A (V k ni)(V k ni) (18) 
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The effect of these two terms is not qualitatively different from the terms K3 and K4 of the main Eq.(g|) and, since 
they depend only on the symmetric product of strains A T ■ A , no new nematic effects should be expected from 
these particular terms. We did not specifically endeavor to obtain the molecular expressions for the corresponding 
constants and, for all practical purposes, we should adopt the approximate equations for the elastic constants ( |l7j ) of 
the general non-uniform nematic rubber-elastic free energy. 



IV. DISCUSSION 



From first glance, the non-uniform elastic energy of a nematic rubber attracts attention by its negatively determined 
structure. One could, perhaps, intuitively accept this qualitative form, — k A 2 (Vn) 2 , as demanding an elastic distortion 
of the sample, subjected to a curvature deformation of the director. One, however, should be careful with conclusions 
because the elastic free energy ( |l5|) contains several matrix products and its overall scalar magnitude depends on 
the mutual orientation of several independent objects: the initial director n° before deformation, the elastic strain 
tensor, the local director n and its gradient. In addition, the effect of AF e i must be considered together with two 
other relevant contributions, the "uniform" nematic rubber elasticity Eq.(|^), which penalizes any deviations of 
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the director with respect to initial n°, and the underlying, conventional Frank nematic elasticity |23|], which is the 
response to any non-uniform deformation Vn. In order to examine the implications of this new elastic energy we 
should consider some simple particular cases. 

In the simplest imaginable situation one clamps the sample, thus prohibiting all uniform deformations. (There 
would still be a possibility for the soft material to deform in a non- uniform fashion with, say, J[A(z) — l]dz = 0, but 
such an effect could be neglected in the first approximation as being of higher order in small deformations) . Taking 
Ay = Sij we arrive at the following expression 

AF e / = + k 3 [1 + (A - l)(n° • n) 2 ] + (A + 2)«a] [(n • V)n] 2 (19) 

-5 [«2 + Ki[l + (A - l)(n° • n) 2 ] + (A + 2)k 7 ] (V fc nO(V fe nj) 
-^«i(A - l)(n° • (n • V)n) 2 - \ 2 {A - l)(n° ■ V fe n)(n° • V fe n) 

+k 5 (A - l)(n° • n)[n° ■ (n • Vn fc )V fc n + n° • (n • Vn)djun] (20) 

Here, in contrast to simple nematics, the difference between n° and n is a director rotation that is penalised even for 
uniform distortions. The energetic cost is substantial, of the order of \x per unit volume and it is thus very likely that 
such a rotation cannot be easily achieved, for example the influence of external electric or magnetic fields is shown 
to be quite insufficient to create a significant deviation of n from n°. If we then take n = n° + <5n(r) and retain, 
as usual, only the leading terms in small non-uniform deviations 8n, the elastic energy simplifies dramatically to: 

AF e i = -~(«2 + k 4 A+ [A + 2]k 7 ) (ion) 2 - ~(k 2 + k 4 A + (A + 2)k 7 ) (n • V x n) 2 (21) 
-^(wi +k 2 + (« 3 + k 4 )A + [A + 2](«e + k 7 )) [(n x V x n] 2 , 

which has an exact form of the Frank elasticity of regular nematic liquid crystals. This is an expected conclusion, 
since we have eliminated all independent vector and tensor variables from our simplified, clamped system leaving 
only gradients of n. What can be considered unexpected is that the remaining combinations of elastic constants yield 
negative square gradient contributions to the free energy! 

It is important to compare the magnitude of the new elastic constants with the standard Frank constants K\\ , K22 
and if 33, which work against the above destabilizing effect. The order of magnitude of our new constants is sufficiently 
easy to estimate because thay are determined by the usual entropic effects of polymers in rubber elasticity. Within 
orders of magnitude we have k ~ pksTa 2 or, in dense thermotropic side-chain systems, k ~ ksT/a ~ 10 J/m. 
This is comparable to values of typical Frank constants. If one considers a main-chain polymer liquid crystal, with an 
extremely large backbone anisotropy ratio A = £i\/£± 3> 1, forming a rubbery network, this comparison can become 
even more favourable to the new negative non-uniform nematic rubber elasticity. 

As we have discussed above, this new elastic energy should be compared with the "uniform" nematic rubber 
elasticity, Eq.(0), which produces for the same clamped case 

F^ = \N x k B T { -^^{5n)\ (22) 

i.e. the magnitude of the director rotation, whether uniform or non-uniform, is penalized by the anchoring to the 
clamped network [^l|j24| l . Let us assume that the destabilizing effect of Eq. has won against the Frank elasticity 
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and consider the overall orientational stability of the sample. For a qualitative analysis we take the angular deviation 
5n ~ #o cos qx and write both the uniform and non-uniform parts of the elastic energy, pulling out common factors: 

F d ~ \N x k B T6l [1 - Na 2 q 2 ] (23) 

(we have neglected all other factors, like A, which we assume to contribute factors of order unity in a side-chain 
polymer system). Within Eq.(p3[), clearly, the orientational instability may occur only with very short wave lengths, 
q^ 1 ~ aN 1 / 2 , which is of the order of network mesh size ({R 2 ) ~ a 2 N). A more detailed analysis of higher order 
contributions, i.e. (Vn) , (Vn) 6 , etc., would then be required. If such small textures were to be formed (in a typical 
experimental case aN 1 / 2 ~ 30 — 50 A), they would be impossible to detect by any optical method and the nematic 
elastomer would appear uniform. It is likely therefore that in the equilibrium state of a nematic elastomer spontaneous 
long wave-length deformations, which would be favoured by the new elastic energy Ai^j, are prohibited by the strong 
uniform anchoring of the director to the network. 

This situation may be changed in three circumstances. 

(i) In a mechanically unconstrained sample, or one with non-uniform distortions such that the sample is globally 
undistorted, there exists a possibility of having so-called 'soft elastic deformations' These do not give rise to 
the uniform part of elastic energy F^, , and happen when a special class of elastic strains is combined with certain 
director rotations in such a way that the network polymer strands are not forced to change their equilibrium shape. 
These strains therefore do not cause a drop in configurational entropy. One example of such a soft strain tensor could 
be A = £ 1 / 2 £ *' 2 (see j2(| for details). Insertion of F$ of Eq.(||), that is Eq.(|J) with uniform fields, shows trivially 
that this class of strains leaves the "uniform" elastic energy unchanged from the unstrained state. It can still allow a 
negative non uniform contribution ~ -K(Vn) 2 . It is expected that networks crosslinked in isotropic state and then 
cooled down into the nematic phase should exhibit such 'soft elasticity' due to their rotational invariance j27j. The 
destabilising negative gradient terms, unrestrained by the positive penalty for uniform rotation, may explain why all 
such elastomers in practice form scattering polydomain textures in thermodynamic equilibrium. 

(ii) Another, perhaps more relevant case, when the destabilising effect of the new non-uniform elastic energy AF e i 
can be felt is during the mechanically driven orientational transitions, for example the ones described in |ll| , pl|] . An 
imposed elastic strain overcomes the barrier to director rotation, given by the uniform energy F^ . The transitions 
predicted and observed were to a uniform, rotated state. But we now see that orientational modulations could occur in 
the strained state. This, of course should take place only when the new constants k are actually more relevant than the 
(stabilizing) Frank constants. We, therefore, predict that materials with higher backbone anisotropy A = £n/£± should 
be more likely to exhibit such spontaneous breaking into orientational domains during various director transitions. 

(iii) It is possible, when there is sufficiently strong surface anchoring of the director in a direction in conflict with the 
principal axes of the imposed strain, that one can have a mechanical Freedericks transition to a non uniform state. 
Even if the new, negative constants do not outweigh the Frank constants in the limit Eq.(|2l|), it is possible that for 
larger A and for n not close to n°, these effects help the nematic rubber elastic contributions overcome the Frank 
penalty and tip the balance in favour of the distorted state with non-uniform nematic textures. 

We appreciate many stimulating discussions with P.D. Olmsted and H. Finkelmann. This research has been sup- 
ported by EPSRC (GV) and by Unilever, PLC (EMT and MW). 
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